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Abstract 

We show that in supersymmetric unified theories such as SO{W), imple- 
mentation of the see-saw mechanism for neutrino masses introduces a new set 
of color triplet fields and thereby a new source of c? = 5 proton decay operators. 
For neutrino masses in a plausible range, these operators are found to have the 
right strength to yield observable, but not yet excluded, proton decay rates. 
The flavor structure of the new operators is distinctive. Proton decay modes 
into a charged lepton, such as £~^7r^, i~^K^ and i^r] where ^ = e or /x, can 
become prominent, even for low or moderate values of tan/? ;S 10, along with 
the VK^ and z/tt^ modes. A distinctive feature is the charged lepton modes 
involving an and/or a /i^ with the ratio T(£'^K^) : T{£~^tt^) ~ 2 : 1. 
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1. Proton decay, if discovered, will constitute impressive evidence for the place- 
ment of quarks and leptons in common multiplets and for the unification of the 
separate gauge interactions of the Standard Model Already, the fact that the 
three gauge couphngs meet at a common scale Mx ~ 2 x 10^^ GeV 0, provided 
they are extrapolated from their measured values in the context of supersymmetry, 
supports the idea of supersymmetric unification. 

Supersymmetric unified theories (GUTs), however, bring two new features to pro- 
ton decay: (i) First, by raising Mx to a higher value as above, they strongly suppress 
the gauge-boson-mediated d = Q proton decay operators, for which e^n^ would have 
been the dominant mode. (In the most straightforward interpretation one obtains 
t{p — > e'^vr") 1(^=6 — lO'^^^^'^ yr., where the uncertainty reflects those from the hadronic 
matrix element and from the masses of the relevant gauge bosons.) (ii) Second, they 
generate d = 5 proton decay operators |0] of the form QiQjQ^Li/M in the superpo- 
tential, through the exchange of color triplet Higgsinos, which are the GUT partners 
of the standard Higgs(ino) doublets. Assuming that the color triplets acquire heavy 
GUT-scale masses, while the doublets remain light, these "standard" d = 5 opera- 
tors, suppressed by just one power of the heavy mass and the small Yukawa couplings, 
provide the dominant mechanism for proton decay in supersymmetric GUT, with a 
lifetime Tp ~ (10^° — 10^^) yr. |^, ^, |^. This range is consistent with present limits, 
and might be within reach of SuperKamiokande. 

The flavor structure of the standard d = 5 operators are constrained by three 
factors: (a) Bose symmetry of the superfields in QQQL/M, (b) color antisymmetry, 
and especially (c) the hierarchical Yukawa couphngs of the standard Higgs doublets. 
Because of these, it turns out that these operators lead to a strong preference for the 
decay of the proton into channels involving a V rather than e"*" or (even) /i+ and those 
involving an s rather than a . Thus they lead to dominant T7K~^ and comparable 
Vtt~^ modes and in some circumstances (i.e., for large tan/5 > 40) to prominent fi^K^ 
mode; but in all cases to highly suppressed e~^7i^ and e~^K^ decay modes. For example, 
in minimal SU{5), with contributions only from the standard d = 5 operators, one 
finds that for tan/3 < 10 (see eg., Ref. i, |, 0): 
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Here i?e7r — Rfnr — 1-2 and Rex — RfiK — 0.12 are the products of the ma- 
trix element and the phase space factors for e^7r° mode etc., relative to the V^K^ 
mode. The factor Ytx refers to the third family contribution relative to the second, 
\YtK\ — \{mtVtdVts) / {'rncVcsVcd)\, and we have neglected any flavor dependence in the 
squark/slepton masses in writing Eq. (1). 

The purpose of this note is to point out that there exists a new set of color 
triplets and thereby plausibly a new source of = 5 operators, in supersymmetric 
unified models like 5*0(10) [Q, which assign heavy Majorana masses to the right- 
handed neutrinos to generate light neutrino masses via the see-saw mechanism 0. 
With a desirable pattern of the neutrino masses, in accord with the MSW solution for 
the solar neutrino puzzle and serving as the hot component of dark matter, these 
new (i = 5 operators are found to compete favorably with the standard ones described 
above. At the same time, the flavor structure of the new operators, related to the 
neutrino masses, appear to be rather universal, and different from the standard ones. 
These new operators allow in general, prominent or even dominant charged lepton 
decay modes of the proton- i.e., p — > i~^K^, and i~^T], where = or /i"*", even 
for low values of tan/3 < 10, along with the neutrino modes p — > T7K~^ and z77r+. A 
distinguishing test of the new mechanism is provided by the prominence of the charged 
lepton modes involving an e"*" and/or a /i"*", together with the prediction r{£~^K^) : 
r(£+7r°) ^2:1. This, as we will discuss, can distinguish the new contributions not 
only from the standard d = 5 operators, but also from certain gauge boson mediated 
effects. 

2. Using standard notations for quark and lepton doublets and also singlets, 
the Yukawa couplings of a color triplet {He) and antitriplet {H'^,) are given by the 
superpotential 
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where i,j are family indices. In the minimal SU{5) model, F and G are the usual 
Yukawa coupling matrices of the standard Higgs doublets with the up and the down 
quarks respectively. To allow for a different flavor structure in the couplings of new 
color triplets, we shall keep F and G general. 

After integrating out He and H'^ superfields, the effective AB ^ superpotential 

is 
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where Mc is the mass of the superheavy color triplet Higgsino {W D McHcH'q) 
The S'f/(3) and SU{2) contractions in Eq. (3) are as follows: 

{QiQj){QkLi) = ea(3-y{ufd'^ - dy'^){ul£i - dlui) 



= e^p,{uri,){ufdr) . (4) 

In terms of component fields, Eq. (3) corresponds to a vertex with 2 fermions and 2 
scalars. For proton decay, the two scalars (which are heavier than the proton) should 
be converted to ordinary fermions by dressing the vertex with a wino or a gluino. 
The contributions of the gluino, which conserves flavor, turn out to be suppressed, 
compared to those of the wino 0], except for the case of large tan/? > 40 (see e.g., 
Ref. 0). For the case of dominant wino contributions, which is what we will mostly 
consider, only the first term {QQ){QL) in Eq. (3) is relevant. 

3. Neutrino masses and new dimension— 5 proton decay operators: Now 
let us identify new candidates for He-, H'q related to neutrino masses. Majorana 
masses for the right-handed neutrinos, which are needed to implement the see-saw 
mechanism, can arise in 50(10) by utilizing the vacuum expectation value (VEV) of 
either a 126// or a 16//. In the case of 126//, one can use the renormalizable coupling 
to matter multiplets (16i) of the form /jj(16il6j)126//; while for the case of 16//, 
one needs to use the effective higher dimensional operator /y(16jl6j)(16//16//)/M. 
Either will generate d = b proton decay operators. We will now discuss each, in turn. 

The case o/126//: In this case, the relevant standard model singlet in the 126// 
that acquires a VEV has the quantum numbers of a di-neutrino "z//jz//j" . This breaks 
S'O (10) to S'[/(5), and as is well known, it has the advantage that it changes {B — L) 



by two units and thereby automatically conserves i?-parity |T0|- Such a symmetry 



neatly forbids potentially dangerous d = A proton decay operators. 

With the 126// acquiring a VEV, there must exist a conjugate 126//, also ac- 
quiring VEV, to cancel the D term. The 126// however has no coupling to the 16, 
because of 5*0(10). The only relevant coupling is therefore 



1^126 = /i, (1646,) 126// . (5) 

Here i, j = 1, 2, 3 refer to generation indices in the gauge-basis. 

The Yukawa couplings fij may be determined (approximately) as follows. The 
light neutrino masses are given by the see-saw formula: m{v^) ~ m(z//^)^/Mj/j, 
where m{v^) denotes the Dirac mass of the ith neutrino, and Mm are related to 
(but are not equal to) the physical Majorana masses of z/j/j. Mj/j are given in 
terms of the matrix elements Mij = /jj ( 126//) as:Q M^/j ^ {Mn, (M11M22 — 



^This pattern emerges if one assumes a hierarchical structure for the Dirac masses without any 
significant hierarchy in the Majorana elements My. 
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Ml^)/Mii, Det{M)/{MuM22-Mf^)}. The successful 50(10) mass relation mb{Mx) = 
mr{Mx) suggests that at least the third family fermions receive their masses primar- 
ily through the Yukawa couplings 16jl6jl0//. This in turn implies that ~ 
mt{Mx) ~ (100 - 120) GeV. The empirical relation m^{Mx) ~ 3m,(Mx) how- 
ever, suggests that dominant contributions to the masses of the second family comes 
from the Higgs component transforming as (2, 2, 15) of SU{2) l x SU{2)r x SU (4)c = 
G224, which contributes in the proportion (1, 1, 1, —3) to the four colors. Such a 
Higgs component, with a VEV of the electroweak scale, may arise effectively either 
from the same 12Qh which gives Majorana masses to the z/)j (see Eq. (5)), or alter- 
natively, and in fact preferably, through an effective operator IGjlGjlO/^ (45//) /M. 
Now 10// X 45// contains the desired submultiplet (2, 2, 15) C 120, which contributes 
only to the off-diagonal mixing (with j = 2, 3), as well as a (2, 2, 1) component. Tak- 
ing both these contributions including the see-saw off-diagonal mixing into account, 
it can be verified that reasonable fits to the second family masses and Vcb can lead to 
~ (3 — 12) X mc{Mx) ~ (1 — 4) GeV. Although not essential for our arguments, 
guided by the masses of u, d, and e, it seems reasonable to take ~ (1 — 10) MeV. 

Thus, with the values of ~ {(1-10) MeV, (1-4) GeV and (100-120) GeV} 
for i = e, fi, T, which are motivated by the observed pattern of masses of the quarks 
and the leptons, one gets, via the see-saw formula:^ rn^uj^) ~ {(^ — 30) x 10^^, (^ — 
5) X 10"'^ and (2 — 3)} eV, provided Mm are nearly flavor universal, within a factor 
of 2 to 3, with MiR ~ M2R ~ M3R ~ (1 - 3) x 10^^ Q^y jg interesting that this 
pattern of masses for the light neutrinos is precisely the one that goes well with the 
MSW solution for the solar neutrino puzzle, involving (z/e — z/^) oscillations (which 
requires m^^ ~ (2 — 4) x 10~^ eV), and with z/,- serving as the hot component of dark 
matter. 

Thus we see that considerations based on quark-lepton masses as well as neutrino 
masses suggest - although they do not mandate - a non-hierarchical pattern for 
the Yukawa couplings of the 126//, with a rather universal Majorana mass Mj/j ~ 
(1 — 3) X 10^^ GeV. This contrasts with the large hierarchy exhibited in the Yukawa 
couplings of the 10// to the three families. 

In the absence of other information, it is reasonable to take the VEVs of all 



relevant Higgs fields (e.g., 126//, 54// and 45//) which break 5*0(10) to the standard 
model symmetry to be nearly equal to the GUT scale, Mx ~ 2 x 10^^ GeV. This also 
ensures that the simple meeting of the gauge couplings within the MSSM framework 



is preserved. With (126//) ~ Mx and M^r ^ /i^ (126//) ^ (1 - 3) x lO^^ Q^y^ 



^These values of light neutrino masses include a reduction of about 50% owing to their running 
from the GUT to the electroweak scale. 
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Figure 1: Superfield diagram inducing the QQQL/M operator in the 126// option. 

get fij ~ (| — 1) X 10^'^. This leads to a strength for the new d = 'b operators (see 
below) of order p jMx ~ (10^^ to 10^^)/Mx, which is of just the right order to yield 
proton decay rate in an observable range. 

To see the origin of the new dimension-5 proton decay operators, let us now 
examine the decomposition of 126 under the subgroup G224 = SU{2)l x SU{2)ji x 
SU{A)c. T26 = (l,3,T0) + (3, 1,10) + (2,2,15) + (1, 1,6). The real (1,1,6) component 
contains the color triplet which we call He and an anti-triplet {H'q). Note that Eq. 
(5) contains the di-quark and lepto-quark couplings of the He and H'f^ C 126// 
respectively (compare with Eq. (2)). Observe that for this case F^j = Gij = fij = fji. 
Thus if a (1, 1, 6). (1, 1, 6) mass term is present, there will be dimension-5 proton decay 
arising from the diagram shown in Fig. |l[ In order to break the SU{5) symmetry 
that is preserved by the VEV of 126//, there must exist other Higgs representations. 
A 45// + 54// is a simple choice. The invariant couplings 

W D A(126//126//54//) + A (T26//T26//54//) (6) 

are then allowed. The 54// acquires a VEV along the (1, 1, 1) direction under G224- 
This supplies the required (1, 1, 6).(1, 1, 6) mass term to induce the new effective 
dimension-5 operator of Fig. Since there is also a 126// 126// mass term, the two 
(1, 1, 6) multiplets coming from 126// and 126// will now mix with an angle parameter 
e. This means that the effective Mq in Eq. (3) is [cos'^e/Mi + 8111^6/1^2] where Mi,2 
are the two mass eigenvalues of the color triplet system arising from 126// and 126//. 

It will be seen later (Section 5) that the interaction of the 54// in Eq. (6) is 
desirable in connection with an attractive mechanism for doublet -triplet splitting, to 
obtain masses for potential Nambu-Goldstone multiplets. Thus there is an intricate 
link between the neutrino masses, doublet-triplet splitting and the proton decay 
operators in this case. 

Only the (1, 1, 6) component of the 126// contributes to dimension-5 proton decay 
operator. It is easy to verify that although the (1, 3, 10) and (3, 1, 10)-components of 
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Figure 2: Dimension 5 proton decay operator involving one neutrino mass vertex and 
one standard color triplet vertex. 



126h contain color triplets, they do not mediate proton decay. This is a consequence 
of SO{10), which prevents the Yukawa couplings of the 126^ to the 16j. 

The case of 16h- In this case, the relevant standard model singlet that acquires 
a VEV has the quantum numbers of a neutrino "z/r". This too breaks 5*0(10) to 
SU{5), but it changes {B — L) by one unit. One can define a Z2 discrete symmetry 
under which 16j's arc odd, but all other multiplets (including 16^^ and 10^^) are even. 
This symmetry can serve as /?-parity. 

Now, with IQh acquiring a VEV, one needs a 16^, acquiring the same VEV, to 
cancel the D-term. S'(9(10)-invariance, together with the Z2 symmetry (under which 
16if is even), now allows the superpotential terms 

W^i6 = ^ftj (16il6,)„ {^hTQh)^ + ^-Qij (16,16,) (16^^16^) + M,,1QhTQh ■ (7) 

Here a in the first term refers to the two possible 5'O(10)-contractions. 

While the first two non-renormalizable terms in Eq. (7) might be taken as quasi- 
fundamental, to be cut off by gravity or string effects at short distances, it is interest- 
ing to examine their possible origins through renormalizable operators. A simple way 
to generate the first term in Eq. (7) that induces the Majorana masses of the vir, is 
via the couplings 16j45j l6if + Mij45i45j. With this couphng alone, which appears 
to be almost inevitable to induce neutrino masses, there are new contributions to 
d — 5 proton decay. 

The relevant diagrams are shown in Figures 2 and 3. In Fig. 2, one of the vertices 
arises from the effective neutrino mass operator, Eq. (7), while the other vertex is 
the standard operator 16jl6jl0i^ proportional to the down quark mass matrix. The 
coupling IGh^Qh^Oh is allowed by all the symmetries, and is also compatible with 
the doublet-triplet splitting mechanism. This term may be desirable since it modifies 
the relation tan/3 = mt/rrib that often occurs in S'O(IO) models, (tan/5 ~ mt/rrib is 
generally problematic for the standard d — 5 proton decay.) This modification arises 
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Figure 3: Another d = 5 proton decay operator from one neutrino mass vertex and 
one standard color triplet vertex. 

because the Y = —1/2 light Higgs doublet, will now be partly from the 10^ and 
partly from the 16^ (see eg., Ref. |12|)- The up-down symmetry preserved by the 
usual WildjlOn Yukawa couplings will now be broken, resulting in tan/? 7^ rrit/mf,. 
Fig. 3 makes use of the fact that the light fermions are not entirely in the 16,, but 
are also contained in 45j, due to the mixing from 16j45jl6/f. 

Analogous to the case of (126h'^, neutrino masses would determine fij (approx- 
imately) by noting that MiR ^ fij (TGhY /M ^ {I - 3) x GeV. Thus we find 
fij (16h) /M ~ (i — 1) X 10"^, with no strong hierarchy in its elements. Conse- 
quently the new contributions from these terms allow charged lepton decay modes of 
the proton to become prominent (see below). 

The IGjlGjlGfflG/f term in Eq. (7) is not directly related to neutrino masses. 
However, it is similar in form to the first term of Eq. (7). It has been used in the 
past to induce realistic fermion masses and mixings in 5*0(10) ]T^. Indeed, note that 
with a single IOh coupling to 16jl6j, the up and down matrices are proportional 
and the CKM matrix reduces to the identity. To correct these bad 5*0(10) relations, 
one needs some additional contributions to the mass matrices. A simple solution 
is to induce the effective operator 16il6jl6i^l6// through the (renormalizable) cou- 
plings IGjlOjlG// + MjjlOjlOj involving superheavy lOj. Then if the mixing term 
16i^l6//10// is also present, so that the SU{2)l doublet from the IQh acquires a 
VEV, the proportionality relations will be corrected and non-zero CKM mixings will 
be induced. From a fit to the masses and the CKM mixing angles, one finds that the 
matrix elements cjij are all of the same order, of the order of the strange quark Yukawa 
coupling, within a factor of 10, gij ~ (10"'^ — lO"'*). These terms would still leave 
the bad relations of minimal 5'[/(5) ms{Mx) = m^{Mx) and md{Mx) = me{Mx)- 
One possible way to correct them, while simultaneously inducing the CKM mixing 
angles, is to utilize the couplings 16jl6 45// + 16^16 lOj/ + IGi/lGj/lOj/ involving 
vector-like pairs of matter multiplets in 16 + 16 (see Fig. 4). This would induce an 
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Figure 4: Possible contribution to 16il6jl6i^l6i^ operator relevant to proton decay 
and fermion masses. 

effective operator (16jl6j)(16j^l6i/)(45ij45ij), which can also serve the purpose of 
the (jij terms in Eq. (7). In this case, again, one sees from a fit to the quark and 
lepton masses and mixings, that the effective couplings gij are nearly universal and 
of order the strange quark Yukawa coupling. 

From all this, it seems natural to assume Qij ~ fij. The d = 5 proton decay 
amplitude (see below) then turns out to have the right strength, as in the case of 
126 H: to yield observable rates. 

Let us focus on the proton decay operators arising from utilizing all three terms 
of Eq. (7). The SU{5) x U{l)x decomposition relevant to Eq. (7) is: 16 = 1"^ + 
5^ + 10~^ where the superscript indicates the U(l)x charges. So the first term in Eq. 
(7) contains the terms 

{k%') (i^iy + (lor^io/O {^h1%) (8) 

while the second term contains 

(5J55) (10-1^^) + (5^07^) (5^1^^) . (9) 

Once the 1% from the 16h and the 1^^ from the 16h acquire (equal) GUT-scale 
VEVs, the first term in Eq. (8) will induce superheavy Majorana masses for the 
zz/j. The second term in Eq. (8) along with the second term in Eq. (9) will lead to 
dimension-5 proton decay. (Recall that the 5^^ + 5^ do not belong to the Nambu- 
Goldstone supermultiplet in the SO{10)/ SU{5) coset space.) The effective operator 
is given by Eq. (3), with the identification of F and G matrices with the / (1%) /M 
and g (ij^^/ /M matrices of Eq. (7) respectively. (The first term in Eq. (9) docs not 



lead to proton decay.) As in the case of 126^, the coupling matrices F and G are 
now non-hierarchical. This turns out to be very significant for proton decay. 

4. Proton decay Modes: We can write down the effective AB ^ four-fermion 
interactions that arise after dressing of the operator in Eq. (3). We do this in 



the mass eigenbasis of the quarks and the leptons. First define the gauge interaction 
eigenstates for the down quarks d[ = Vijdj and for the neutrinos u'- = Vj'uj with V 
and being the CKM matrix in the quark sector and the lepton sector. Then the 
effective AB ^ four-fermion interaction is 

X [«df)(4M)(/«,«.) + /K,^r 

+ «"«f)(rf7z/;)(/(«„4) + /K-,£r; 
+ «rff)(«j£r) + ^D)] . (10) 

Here the color indices (a,/5,7) and the fiavor indices {i,j,k,l) are understood to 
be summed over, and the fermion fields paired together in parentheses are spin- 
contracted to singlets. / is a loop integral, with magnitude of M^gy, defined as 

/(a, b) = ( -[a-, b]] . (11) 

Note that all the matter fields in Eq. (10) belong to weak isodoublets. Traditionally 
it has been argued that the third term dominates; one of our main points here has 
been to emphasize that this need not be so. 

Since Eq. (10) is in the mass eigenbasis, the Yukawa couplings Fij and Gij are 
not the same as Fij and Gij of Eq. (2). F remains symmetric, it is related to F by 
F = VJFVu while G = (VjGVujV , where Vu is the unitary matrix that rotates the 
left-handed up quarks in going to the mass eigenbasis. V is another unitary matrix 
that parameterizes the mismatch between the up quark and the charged lepton mass 
matrices |13|, V = V^Ve. Note that G is not symmetric. 

The proton decay rate and branching ratios can be obtained from Eq. (10). As 
we have emphasized, several considerations suggest the matrices F and G may not be 
hierarchical, so proton decay into alternative fiavor modes could a priori have similar 
rates. Consider first the decay into charged leptons. The relevant interactions are the 
second and the fourth terms in Eq. (10). In the second term, we must put j = k = 1 
since the operator has to have only u quarks. Flavor antisymmetry then requires 
that z = 2 or 3. Similarly in the fourth term, i = j = l,k = 2,3. Noting that 
d'2 ~ Vcdd + VcsS and d'^ ~ Vtdd + VtsS for proton decay, the amplitude for charged 
lepton decay can be written as 



A{p ^ ^+) (X {{u"d''){u^er) [V,d {GiAi - G2iFu) + Vu {GuF^i - 

+(M"/)(n^£r) [Vcs {GuF2i - GszAi) + V^ {GuF,, - Gs^Ai)] } • (12) 
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Now if the matrices F and G have no strong hierarchy in their flavor-dependence, 
then the terms proportional to Vtd and Vts can be ignored in Eq. (12). The error 
introduced is only of order ~ 1/20, (A^ = smOc — 0.22). This observation leads 
to predictions for the branching ratios of certain (in general) prominent modes, which 
are independent of the flavor structure in F, G\ 

r(p ^ i+K^) 




fiin'Oc 1 - —I R- (13) 



Here = or and R is the ratio of the two relevant hadronic matrix element- 
squared. The chiral Lagrangian estimate for R given in Ref. ||^ is 



R 



[l + D + F)\ 



2 



(14) 



Using D = 0.81, -F = 0.44 for the chiral Lagrangian factors and with tub = ttia = 
1150 MeV (as in H), we obtain ~ 5. Thus 

r{p £+7r°)/r(p ^ i+K^) ^ 0.5 . (15) 

Another interesting mode is p ^ i'^r] for which one has 



r(p-.£+7rO) V ml) \{l + D + F)[ 



0.35 . (16) 



While all charged lepton modes are expected to have similar rates, the ratios 
such as T{p —>■ e+7r'^)/r(p —>■ fi'^n^) cannot be predicted quantitatively since they are 
sensitive to the flavor structure of F and G: 



T{p ^ e+7r°) 
T{p /i+vrO) 



{G11F21 — G21F1 



{G12F21 — G22FU 



2 



(17) 



If one uses 16h to generate the Majorana masses, the matrices G and F are 
independent, so the ratio in Eq. (17) is in general expected to be of order unity. 
However, if a single 126 h is used for this purpose, one has the asymptotic relation 
G = FV. The flavor-dependent renormalization of this relation is small. If in 
addition the off-diagonal entries in the Jarlskog matrix V are small, then one has 
a cancellation in the amplitude for e~^iT^. In this case, one would expect the n'^n^ 
mode to dominate over the e+vr^ mode. Note that in the 16h option, however, since 
G and F are independent, the two modes are expected to be comparable even if V 
has small off-diagonal entries. 
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In the case where 16h is used to generate the see-saw neutrino masses, as noted 
earher, new d = 5 operators arise simply from the first term in Eq. (7). Then the 
factors Fij, being related to the masses, are non-hierarchical, while the factors 




dropped. Again one finds that r(//+X°) : r(/i+7r°) ~ 2 : 1, as in Eq. (15). A 
similar remark apphes to r(e+ir°) : r(e+7r°). Now, however, A{e'^K^) oc (F21 + 
sm9cFii)md, so that r(/i''"K°) is expected to be considerably larger than r{e'^K^) 
(barring fortuitous cancellations). 

It is worth noting that, independent of the relative importance of positron modes, 
the /i+TT^ and fi~^K^ modes arising through the new d = 5 operators can still compete 
favorably with or supersede the VK^ modes, for even small or moderate values of 
tan/5 < 10. By contrast, for the standard d = 5 operator, the fi^n^ and fi^K^ modes 
can be prominent only for very large tan/5 > 40 through gluino dressing. Thus a 
study of these decay modes and determination of tan/3 could distinguish between the 
standard and the new d — 5 operators, even for the case of 126h- 

The decays p VK'^ and p — > Vn^ from the new dimension-5 operators have rates 
comparable with the charged lepton modes. The expectation for V{p — > VK^)/T{p — > 
Vti^) is similar to but not exactly the same as the case of minimal supersymmetric 
SU{5). The flavor structure relevant for the V decay is given by (in the limit of 
neglecting terms of order < ~ 1/20): 



A{p ^ 17,) oc {(«"d'^)(sV,)[KdK.(AiG2£ - FuGu) - VMF22GU - F,2G2l)] 

- {u''s''){d''ue)[V,dVUF22Gu - ^2^2^) - ^.^^(^162^ - F^Gu] 

+ {u"d''){d^ue)[V^,VUFnG2e - ^2^1,) - ^^(^22^1, - F,2G2e)]}. (18) 



Note that in Eq. (18), the VudVcs term (the first term) has no mixing angle suppres- 
sion, but all the remaining terms are suppressed by factors of at least A^. It turns 
out that the matrix element for the {u'^s^){d'^vi) is suppressed by a factor ~ A^, 
compared to the {u'^d^){s"'vi) term, leading to an overall suppression factor of ~ A^ 
in the terms in the second line. Neglecting these sublcading terms one finds 



This prediction holds if the contributions from the new d — b operators dominate. The 
analogous number for the standard d — h operators in minimal SU{b) is near |, but 
there is considerable uncertainty in this case owing to possible cancellation between 




Gij ~ V*jm'^j exhibit a hierarchy. The charged lepton modes become prominent even 
in this case. For the amplitudes, one obtains: A{^'^K^) oc (Fn — sin^c'-^2i)"T's; while 
A{iJi^T^^) oc sin^c(Fii — sin^c-^2i)"^s, where terms of order A^ ~ 1/20 have been 




(19) 
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the second and the third generation contributions. It is difficult to make quantitative 
estimate of ratios such as T{i'^K^)/T{vK^) etc., arising from the new d = 5 operators, 
since there is some flavor dependence. Comparing the flavor structure in Eq. (12) and 
Eq. (18), one infers that the rates for both these modes are similar, with the VK^ 
mode slightly preferred over the l^K*^ mode, owing to a matrix element enhancement 
and the availability of final states with all three neutrino flavors. However, there are 
terms that differ in the two amplitudes which are of order Ac — 1/4.5, and so it is 
difficult to make a more precise quantitative estimate. 

The decay rate of the neutron can be obtained from a general operator analysis 
g. For example, T{p £+7r°)/r(r2 f+vr") ~ \ and T{n VK°)/T{p VK+) ~ 



1.8. 

5. Doublet— triplet splitting and the standard = 5 proton decay: An 

important issue that any realistic GUT model faces is the question of doublet-triplet 
splitting. While the SU{2)l doublets in the IQh of SO{10) have to be light in order 
to trigger electroweak symmetry breaking, their color triplet GUT partners have to 
remain heavy at the GUT scale, since they mediate proton decay. One attractive fea- 
ture of supersymmetric S'O(IO) is the existence of a natural doublet-triplet splitting 



mechanism [T^, 16]. This mechanism can bring in a numerical suppression in proton 
decay arising from the usual (16il6j)10H coupling that gives rise to the quark and 
lepton masses. These operators, which lead to p — » VK^ as the dominant mode, are 
somewhat problematic in supersymmetric SU (5) since the predicted rate is near the 
experimental limit 0, |^. 

The doublet-triplet splitting mechanism in 5*0(10) utilizes the superpotential 
couplings 

Wdt = A10h4:5h10h + M'IO'jjIO'h , (20) 

where the IOh is the Higgs superfield which contains the two Higgs doublets of MSSM 
and where lO'fj is another field with a GUT scale mass M'. Once the 45h acquires 
a VEV along the {B — L) direction, (45^^) = diag.(a, a, a, 0, 0) ® T2, the color triplet 
mass matrix M. and the S'f/(2)-doublet mass matrix JVi' become 

M^{? t): M'^il 1). (21, 



Aa M'J ' VO M' _ 

This gives GUT scale masses to all the triplets from the IOh and 10^ while one pair 
of Higgs doublets from 10^ remains light. Proton decay amplitude mediated by the 
color triplets in 10^ is now proportional to {A4^^)u = (M'/A^a^), so by choosing M' 
somewhat smaller than Aa, one obtains a numerical suppression of proton decay. M' 
cannot be too small compared to the GUT scale, however, since that would result in 
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a large positive contribution to the predicted value of a^lMz). The shift in a^lMz) 
from the doublet-triplet sector alone is given by 

AasiMz) = -^^^^^l\nUM-')^,Mx] • (22) 
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If {M-^)u = [10^^ GeV]-\ Aas = +0.005, which might be acceptable. If {M' 
[10^^ GeV]^^ so that proton decay from this operator is unobservable, then Aas{Mz) — 
0.011 which seems excessive and would require a cancellation from some other thresh- 
old effects. This implies that p — * VK~^ cannot be suppressed to an unobservable level, 
at least in this simple doublet-triplet splitting scheme. 

In order for the VEV of 45h to be along the {B — L) direction to a great ac- 
curacy, the 45 H should not couple or should couple only weakly to the 16 h + IQh 
sector. There is a danger of having pseudo-Goldstone multiplets, which could upset 
the unification of the gauge couplings, if such cross-couplings are prevented. The po- 
tential pseudo-Goldstones are the {(3, 2, i) + (3*, 1, — |) + H.c.} components, under 
SU{3)c X SU {2)l X U{1)y, from the 16h and the 45h- This issue has been addressed 
in the case of using 16h + 16// to break 5*0(10) to SU{5) [|T6|. Here we point out 



that the same mechanism will work if the 16// + 16// is replaced by a 126// + 126//. 
The couplings used in Eq. (6) that were relevant for the proton decay operators are 
precisely the ones that can give masses to all these would-be-Goldstone multiplets. 
So in this case the proton decay amplitude is closely related to the doublet-triplet 
splitting mechanism. 

6. Gauge boson— mediated d = 6 versus the new d = 5 proton decay 
operators: So far we have focussed on the d = 5 proton decay operators. In the 
simplest supersymmetric GUT models, the gauge boson mediated d = 6 operators are 
suppressed relative to the d = 5 operators. However, enhancement of the d = 6 oper- 
ators could occur for a variety of reasons. Possibility arises in flipped SU{5) x f/(l) 
17 1, in non-supersymmetric two-step breaking of 5*0(10) models [jl8|, or possibly 



even in supersymmetric SU{5) [T^ and 5*0(10) models with large threshold correc- 
tions, wherein the relevant (X, Y) gauge boson masses are of order 10^^ GeV. Ex- 
change of these particles would lead to a dominant e^vr^ mode, with proton lifetime 
~ (10^^ — 10^^) yrs., compatible with current limits. We wish to point out that should 
the e^TT^ decay mode of the proton be observed, one can empirically decide whether 
it has its origin in the gauge-boson-mediated c? = 6 or in the new d = 5 operators 
discussed here. For the gauge-mediated case, e^K^ will be strongly suppressed com- 
pared to the e'^vr" mode, by the Cabibbo angle (sin^6'c — 1/20), phase space (~ 1/2) 
and relevant matrix element- squared (~ 1/2), so that [r(e+-ft'°)/r(e^7r°)]rf=6 ~ 1/80. 
By contrast, for the new d = 5 operators, we have shown that e~^K^ exceeds e'^ir^ 
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rate by about a factor of two. 

7. In conclusion the following remarks are in order. 

(i) While the new d = 5 proton decay operators seem to be best motivated by 
their link to neutrino masses, we wish to note that the results presented here are 
more general. Indeed, the prominence of charged lepton modes and the predictions 
for certain branching ratios depend only on the assumed non-hierarchical nature of 
F and G in Eq. (10). This condition might be satisfied in other contexts as well. For 
example, if the CKM angles are induced in 5*0(10) by coupling the fermions to two 
IOh of Higgses, a reasonable fit may be obtained when the SU{2)l doublets and the 
color triplets from one of the lOjy has non-hierarchical couplings, with strength of 
order the strange quark Yukawa coupling. The exchange of these color triplets would 
lead to results similar to the ones presented here. 

(ii) As has been discussed in the literature, ci = 5 proton decay operators such as 
QQQL/M, belonging to 16il6jl6fel6i/M, could be induced not only by the exchange 
of GUT-related color triplets, but also by other effects including exchange of the heavy 
tower of color triplet string states. These are allowed by 50(10) symmetry as well 
as as i?-parity (or the Z2 symmetry mentioned in Sec. 3). In any supersymmetric 
theory, these non-renormalizable operators must somehow be suppressed at least by 
a factor of 10"'' (if M ~ Mpi) in order not to conflict with observed limits on proton 
lifetime. For a discussion of this issue and its possible resolutions through the use of 



flavor symmetries, in the context of string-derived solutions, see Ref. |2^, and in a 
non-string context, see for example, Ref. pi . 



(iii) It is worth noting that there may be circumstances where the potentially dan- 
gerous GUT-related color triplets are projected out, e.g., below the compactification 
scale of a string theory that leads to a non-GUT symmetry like G224 C 5*0(10) 
but the components of \Qh and IQh providing Majorana masses of the right-handed 
neutrinos may still exist below the string scale. In this case, the standard ci = 5 
operators will be absent, but the new d = 5 operators discussed here could still be 
effective. 

In summary, we have shown that in a class of grand unified theories including 
supersymmetric 50(10), there can be a significant link between the neutrino masses 
and proton decay. In the process of generating neutrino masses one typically induces 
a new source of d = 5 proton decay interactions with an interesting strength. The 
flavor structure of these new d = b operators is distinctive. In contrast to the standard 
d = 5 operators, the new ones can lead to prominent (or even dominant) charged 
lepton decay modes, such as l^ir^^i'^K^ and t^rj, where i = e or fi, even for low or 
moderate values of tan/3 ^ 10, along with VK^ and z/tt"*" modes. A distinguishing 
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feature of the new mechanism, relative to d = Q vector exchange, is the predicted 
ratio T{i+K'^) : r(^+7r°) -2:1. 
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